Microrheology, stress fluctuations and active behavior of living cells 
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We report the first measurements of the intrinsic strain fluctuations of living cells using a recently- 
developed tracer correlation technique along with a theoretical framework for interpreting such 
data in heterogeneous media with non-thermal driving. The fluctuations' spatial and temporal 
correlations indicate that the cytoskeleton can be treated as a course-grained continuum with power- 
law rheology, driven by a spatially random stress tensor field. Combined with recent cell rheology 
results, our data imply that intracellular stress fluctuations have a nearly 1/uJ 2 power spectrum, as 
expected for a continuum with a slowly evolving internal prestress. 
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An accurate physical picture of the viscoelasticity and 
motion of the cytoskeleton is crucial for a complete un- 
derstanding of processes such as intracellular transport 
FiJ , cell crawling , and mechano-chemical transduction 
[3 • Microrheology 0] , based on the analysis of embedded 
tracer particle motion, has recently emerged as an exper- 
imental probe of cytoskeleton viscoelasticity and dynam- 
ics 0, IE • The viscoelastic properties of eucaryotic 
cells arise from an intricate network of protein filaments 
driven by specialized motor proteins and directional poly- 
merization, that convert the chemical energy of adenosine 
triphosphate (ATP) to mechanical work and motion. A 
cell is thus a nonequilibrium soft material whose fluctua- 
tions are actively driven. Unlike the thermal fluctuations 
in an equilibrium material, the amplitude and spatial dis- 
tribution of active fluctuations can be controlled via bio- 
chemical signaling pathways; perhaps allowing the cell to 
locally adjust its' mechanical properties to suit its' needs. 
Indeed, microscopic force generators play a central role 
in existing cell mechanics models such as the sol-gel p|, 
soft glassy rheology Q and tensegrity hypotheses. 

In this Letter, we extend a recently introduced method, 
termed two-point microrheology 10], and show that it 
can be used to characterize the activity of intracellular 
force generators by directly measuring a cell's intrinsic, 
random stress fluctuations. Our experimental data and 
theoretical framework show that a cell can be modelled as 
a coarse-grained viscoelastic continuum driven by a spa- 
tially random stress field having a 1 /uj 2 power spectrum 
in our observable frequency range, I < uj < 60 rad/s. 

There are two distinct approaches to microrheology: 
the active approach measures the displacements of tracer 
particles induced by external forces and the passive ap- 
proach measures fluctuations of particle positions in the 
absence of driving forces. The active approach provides 
a direct measure of the complex shear modulus h{uj). 
In equilibrium systems the passive approach also mea- 
sures fi(uj) because of the fluctuation-dissipation theo- 
rem (FDT) Literature results in cells using single- 
particle versions of the two approaches yield shear moduli 



differing by orders of magnitude and exhibiting qualita- 
tively different frequency dependencies 0, ||| . These re- 
sults are further limited by ambiguities associated with 
tracer boundary conditions and medium heterogeneity 
and, more importantly, the applicability of the FDT. 
We show theoretically that data from passive two-point 
microrheology and active response experiments can be 
combined to measure the activity of molecular motors 
- specifically, the power spectrum of intracellular stress 
fluctuations - despite the cell's heterogeneity. We will 
first discuss current microrheology approaches and our 
measurements . 

Microrheology relies on the Langevin equation for 
the velocity v of a tracer particle, 



md t v(t) 



dt'"/(t-t')v(t') + f E (t)+f R (t), (1) 



where m is the mass of the tracer particle, fjs(i) is 
the external force, and ffl(i) is the random force aris- 
ing from the medium. Stokes law states that 7(w) = 
67ra/i(u;)/(— ioj) |12| for spherical tracers with radius a 
and no-slip boundary conditions in an isotropic, homo- 
geneous viscoelastic medium with a complex shear mod- 
ulus fJ,(oj). In the active scheme, fi?(t) can be set to 
zero, and the displacement is r(ui) = x(u>) fe(a;), where 
x(w) = {— ioJ [—imoj + j(oj)]} is the response function. 

Several measurements of x( w ) to extract the elastic 
moduli of the cytoskeleton of living cells have been per- 
formed 0,0. In particular, Fabry et al. have reported 
accurate measurements of the linear response to an ap- 
plied torque of a few-micron-diameter magnetic bead at- 
tached to the outside of a cell and strongly coupled to 
the cytoskeletal network via the cell's integrin receptors. 
Remarkably, they found that over five decades in fre- 
quency, from 10~ 2 to 10 3 Hz, the shear modulus of sev- 
eral cell types was of order 1 kPa with a power law form 
~ uj 13 with 0.10 < /3 < 0.30. 

Passive 'one-point' microrheology is based on the cor- 
relation function C TiT - (t) = (ri(t)rj(0)}. In thermal equi- 
librium, this correlation function is related through the 
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FDT to the response function x(w): 

The final form is valid in the window uj a < to < ujb, with 
the lower frequency u a ~ I0~ 2 Hz set by the compres- 
sional mode of the network and upper frequency u>b ~ 10 5 
Hz by inertial effects |T^j . Experiments typically measure 
tracers' mean-squared-displacement (MSD), (Ar 2 (r)} = 
(An{t,T)Ari(t,T)), where Ar;(t,r) = r;(f + r) - r<(t) 
and the brackets represent time and ensemble aver- 
ages. (Ar 2 (r)) is related to the correlation function by 
(Ar 2 (r))=/f (l r e-^)C ri rM- 

Two-point passive microrheology experiments mea- 
sure a two-particle displacement tensor Dij(R,r) = 

(Ar{ (t, r)Arj 2 (i, r)), where the superscripts identify 
the two different tracers, the brackets represent an en- 
semble average over all tracer pairs and time, and R is 
the separation between the two tracers. In thermal sys- 
tems, Dij(R,r) can be related via the FDT to the shear 
moduli by two-particle equivalents of C riTj (uj), but with 
the l/a scaling replaced by l/R scaling. Importantly, it 
has been shown [l3l that Dij(R,r) is independent of 
tracer size, shape and boundary conditions (to leading 
order in l/R). In cells, this allows the use of endogenous 
particles, unlike existing methods that attach or inject 
synthetic tracers which may perturb the cell. In systems 
where Eq. J2J is valid, Dij(R,r) and (Ar 2 (r)) have the 
same functional form vs. r, differing only by a geometri- 
cal constant. This leads to the definition of a two-point 
MSD, which may be thought of as an ideal Stokes par- 
ticle advected by random fluctuations of the medium: 
(Ar 2 (r)}2 = (2R/a)D rr (R,T), where D rr is the tensor 
component of Dij parallel to R. 

We made passive microrheology measurements on two 
cultured cell lines for which shear moduli had been 
measured previously Q. All cells were cultured in Del- 
beco's Modified Eagle Medium supplemented with 10% 
calf serum and 50 mg/ml gentamicin in 5% CO2. F9 
cells were cultured in gelatin-coated tissue-culture flasks. 
Cells were passed into glass-bottom cell-culture dishes 
with collagen-coated coverslips and allowed to incubate 
for either 4-6 hours (F9) or overnight (J774A.I) prior to 
experiments. For tracers, we visualized endogenous re- 
fractive particles in the cells using differential interference 
contrast (DIC) microscopy and computed one- and two- 
particle correlations from the same multi-particle video 
tracking trajectories [l5|. As our analysis assumes a 
three-dimensional continuum, we excluded trajectories 
from the thin lamellar region of the cell and the mechan- 
ically distinct nucleus, focussing instead on the midplane 
of the 6-8/zm thick cells. We typically observed ~ 100 
sub-micron tracers (presumed by morphology to be lipid 
granules and mitochondria) within our I /jm focal depth. 
Thirty minutes of data yielded ~ 10 7 tracer positions per 
cell, with respective time and space resolution of 1 /60th 
sec and 20 nm. 
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FIG. 1: One (closed symbols) and two-particle (open sym- 
bols) displacements vs. lag time r for endogenous tracers 
in two types of cultured cells, J774A.I mouse macrophage 
and F9 mouse carcinoma. Super-diffusive behavior, exponent 
greater than 1, indicates the effect of non-thermal fluctua- 
tions, (inset) RD rr (R,T = 0.067 sec) in units of 10 _4 /xm 3 , is 
nearly constant in both cell types, as expected for a coarse- 
grained isotropic continuum. 



Typical MSD data for cells are presented in Fig. ^ 
Although visual examination of the images occasionally 
shows ballistic tracer motion (trafficking), driven presum- 
ably by kinesin and dynein motor proteins along micro- 
tubules, the ensemble averaged (Ar 2 (r)} is dominated 
by random, apparently diffusive motion. From the two- 
point data, we first verified that D rr ~ I /R, [Fig. 2] (in- 
set)] in the accessible range 2 < R < 8/xm, comparable 
to results in other studies 0, 0. In thermal systems 
at least, this finding suggests that the medium may be 
treated as a coarse-grained homogeneous continuum on 
the scale R. Strikingly, the two-point displacements of 
intracellular tracers exhibit a super-diffusive behavior: 
(Ar 2 (r)} 2 ~ t q with exponent 1.30 < a < 1.60. If a 
living cell were an equilibrium viscoelastic medium with 



as reported in Ref. Eq. (|2J) would im- 



ply that (Ar 2 (i)} — t 13 . Our data show that a > (3 
and, thus, the FDT breaks down. While the one-point 
MSD suggests that the motion of single endogenous trac- 
ers is dominated by diffusive motion relative to the net- 
work, which may be thermal or non-thermal in origin, 
recent experiments have shown (Ar 2 (r)) ~ r 15 for large 
(several pan) phagocytosed tracers |6j . Such one particle 
measurements, however, can not distinguish between ad- 
vection by a driven continuum or trafficking relative to a 
stationary network. In constrast, our two-point MSD re- 
sults unambiguously indicate that the cytoskeleton itself 
has large strain fluctuations driven by nonthermal forces. 

In living cells, active motors not only modify the vis- 
coelastic response of the cytoplasm, they also give rise to 
random, non-thermal stress fluctuations that cause tracer 
particles to be subjected to random nonthermal forces. 
In such a system, Eq. J2J must be modified. We will out- 
line a careful derivation below, but we can easily guess 
it as follows: In thermal systems, the stress tensor has a 



3 



-13 



10 

io" 14 U 



^ io- 15 L- 



I 10 

< 



-16 



10 



-18 = 



- 1 1 1 — I — I I I I 



- 1 1 1 — I - 



F9 

J774A.1 - 




—I I I I I I L 



I I I 1_ 



10 

co (rad/s) 



100 



FIG. 2: The intracellular stress fluctuation spectrum A(cj) 
as measured in two-point microrheology. The cell lines are 
J774A.1 macrophage and F9 carcinoma. The dashed line rep- 
resents the thermal spectrum, 2fcsT '//' '(lj) / 'to , where is 
measured in Ref. 0|. 



random component Sij (x, t) whose local fluctuations are 
described by 

(s li (x,w)s fci (x',w')) = 2ttA{uj)5{u + lu')5 3 (x - x') 



(3) 



where A(lu) = 2k B Tp" {lu)/uj as dictated by FDT fPf\ . 
In non-equilibrium systems, there are stress fluctuations 
of the same form, but with A(w) not locked to /j,"(oj) and 
we expect 



CviTj ~ Si 



A(w) 



67ra|/i(w)| : 



(4) 



where A(lu) is now interpreted as the power spectrum of 
the stress fluctuations whose microscopic origin is the ac- 
tivity of the motors. Equation (QJ suggests that (i) trac- 
ers can exhibit superdiffusive behavior provided A(uj) di- 
verges sufficiently at small w, and (ii) if an independent 
measure of fJ,(u>) exists, then A(u>) can be sensibly ex- 
tracted from passive correlations data. 

To avoid artifacts associated with medium heterogene- 
ity and non-Stokes' boundary conditions, we used Eq. © 
below, which is the two-point equivalent of Eq. (0J, along 
with rheological data reported in Ref. |4| , to convert our 
two-point data for both cell types to A(uj) as shown in 
Fig. We found a nearly 1/u 2 spectra. These typi- 
cal results were replicated on 8 cells of each type. The 
variations of a and the power spectrum exponent among 
cells of each type, and within different regions of a single 
cell, were comparable within our measurement error. A 
1 /uj 2 spectrum corresponds to a linear decay in time of a 
stress-stress correlation function within our experimental 
time window, and would be a natural consequence of slow 
evolution of intracellular stress [Ts| . Note also that A(w) 
is a few orders of magnitude greater than the thermal 
spectrum (dashed line), affirming that these fluctuations 
are driven by a nonthermal mechanism. 



In the case that the shear modulus and the stress 
spectrum have power-law forms, our result can be sim- 
ply stated. In living cells, with 'power-law' modulus 
~ ojP and a A(oj) ~ u;~ 7 spectrum, Eq. Q im- 
plies that (Ar 2 (r)) ~ T 7+2/3-i^ jf ^he 1/uj 2 spectrum we 
observe is universal, then one can measure the rheology 
exponent (3 via passive two-point measurements, using 
the formula (Ar 2 (r))2 ~ r 1+2/3 . The correspondence of 
our explicitly intracellular results with those of Fabry et 
al. 3 supports their assertion that they are measuring 
intracellular rather than cortical viscoclasticity. 

Now, we turn to our theoretical task and address 
the fundamental question: Given that living cells are 
highly heterogeneous non-equilibrium systems in which 
the FDT does not apply, do two-point passive microrheo- 
logical experiments overcome ambiguities associated with 
tracer boundary conditions and medium heterogeneity, 
and in conjunction with response measurements, extract 
the power spectrum of the continuum stress fluctua- 
tions? We first postulate that a living cell at large length 
scale effectively is an incompressible, viscoelastic medium 
characterized by a displacement field u(x, t), whose equa- 
tion of motion is given by 



P 



d 2 Uj 
dt 2 



dj<Jij{-x.,t) 



(5) 



where p is the coarse-grained mass density, (x, t) is 
the non-thermal random force density arising from the 
motors, and <xy(x, t) is the stress tensor. Within lin- 
ear response theory, the stress cr^ (x, t) is related to the 
strain by Oij (x, to) = 2/z(w)wi 3 (x, w), where Uij = (djUi + 
diUj)/2 is the strain tensor and (j,(cj) the frequency- 
dependent shear modulus. Since there should be no net 
external force in a stationary cell, the random force must 
take the form /j(x, t) — djSij(x,t), where Sij(x,t) is a 
random stress tensor. We assume that its average is zero 
and variance given by Eq. JJjJ . To relate to our microrhe- 
ology experiments, we relate one- and two-particle corre- 
lation functions to stress fluctuations A(u>) and response 
to external forces fj,(ui) as follows. 

We consider two tracer particles of radius a placed in 
this random medium and ask: What forces are exerted 
on each particles? For simplicity, we assume that the 
heterogeneities near the vicinity of the tracers are well 
reflected by a local effective shear modulus fJ,*(u>) which 
may be different from in the bulk. Let particle 1 

(2) at x (x') undergo a displacement s^\lu) (e\ 2 \lu)). 
First, we decompose u(x, t) into an average part and 
a fluctuating part: u(x, t) — u(x, t) + u(x, t), and 
solve Eq. (0 subject to boundary conditions: u»(|x| = 

a, u>) — s\ n \u>) and Uj(|x| = a,ui) = 0, on the surface 
of the nth particle and Uj(|x|,a;), Uj(|x|,a;) — > 0, far 
away from the particles. The total force exerted by the 
medium on each particles has two components: an aver- 
age, F^ n \uj) — J s dS hj (tij(x,u)), and a random parts, 

Pj; n \uj) = J s dS hj 5ij(x, lj), where = 2/J,(u>)uij and 
ijij = 2fi(Lo)iiij are, respectively, the average and fluctu- 
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ating stress, and hj is the unit surface normal pointing 
towards the center of each particle. It is straightforward 

to compute F^ n \ui) = J2 m x7/ ( w ) e i C^O) where 

Xij ( w ) is the two-particle response matrix given by 



t, — \ and \1]' 2 \lo) 



to the lowest order in 1/R, where R = |x — x'|, the 
distance between the two particles and R = R/i?; the 



noise correlators are given by (F^ n \uj) F^" l> {— lo)) = 
A(w) Xij 1 ^ n,m \— tj)/fx(— u>) 19]. These results imply 



that {e[ 1 \uj)ey (—w)) depends in a complicated way on 
H*(to) and In contrast, the cross-correlation func- 

tion as measured by two-point microrhcology is 



-(i) 



D rr ( J R,a;) = (4 1) (-)4 2) (-)) = ^## 



6nR \h{uj)\ 



(6) 



to the lowest order in 1/R, which depends only on /i(uj) 
and A(u>) in the bulk, and is independent of the tracers' 
size, shape or boundary conditions. Thus, apart from a 
geometrical factor, Eqs. and © are equivalent, and 
this justifies our interpretation of our experiments. Fur- 
thermore, the 1/R behavior shown in Fig. ^ (inset) not 
only implies that the living cell can be treated as a con- 
tinuum, but also requires that the random force fi(x,t) 
arises from a stress tensor of the form given in Eq. (J3j). 

Lastly, we propose a microscopic picture of motor ac- 
tivities and derive a fluctuating stress tensor. Since mo- 
tors are small but finite objects, their activities disturb 
the ambient cytoplasm in the form of a point dipole 
pof. Since there is no net external force inside a sta- 
tionary cell, by Newton's 3rd law, the force exerted on 
the cytoplasm by a motor must be equal and opposite to 
that of the cytoplasm on that motor. Thus, the part of 
the stress tensor arising from deviations 5c a (x, t) of the 



coarse-grained activity density of motors from its average 
c a is Sy(x, t) = r (hihj — \8ij) 5c a {x, t), where n is the 
direction of the point dipole and T is the energy scale. 
Since there is no preferred direction, we can average over 
all angles to obtain 



(sij(x,cj)s w (x', -w)) = 



15 



SikSji + SuSjk — - SijSu 



x (Sc a (x,cj)5c a {x',-u)}, (7) 

which has exactly the same form of Eq. J3J if the fluctu- 
ations of motor activity is delta correlated in space. Fur- 
thermore, we can estimate the power associated with the 
stress fluctuations by computing the power dissipated per 
tracer: (P) = J ga u fi" (J) / |^H| 2 ~ KD 16 W. 
Estimating ~ 10 3 tracers per cell, the power to drive one 
cells' fluctuations is of the order of 10 _13 W, which, as 
1% of typical cell metabolism, is not unreasonable. 

In conclusion, the interpretability of two-point mi- 
crorhcology in complex media, including cells, serves 
to clarify an otherwise confusing set of biophysical ob- 
servations, supporting efforts to model cells as three- 
dimensional continua rather than cortical shells and sug- 
gesting that the cytoskeleton is a highly dynamic, ac- 
tively stressed network. Future work to extend the tem- 
poral range and statistical power of such measurements 
should enable mapping and non-trivial spectroscopy of 
intracellular rheology and stress. Finally, such techniques 
may prove fruitful in other systems with non-thermal 
fluctuations, such as granular and jamming media. 
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